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ABSTRACT 

The chromosphere of the quiet Sun is a very important stellar atmospheric region whose thermal 
and magnetic structure we need to decipher for unlocking new discoveries in solar and stellar physics. 
To this end, we need to identify and exploit observables sensitive to weak magnetic fields (B<100 G) 
and to the presence of cool and hot gas in the bulk of the solar chromosphere. Here we report on 
an investigation of the Hanle effect in two semi-empirical models of the quiet solar atmosphere with 
different chromospheric thermal structures. Our study reveals that the linear polarization profiles 
produced by scattering in the Ca ii IR triplet have thermal and magnetic sensitivities potentially of 
great diagnostic value. The linear polarization in the 8498 A line shows a strong sensitivity to inclined 
magnetic fields with strengths between 0.001 and 10 G, while the emergent linear polarization in the 
8542 A and 8662 A lines is mainly sensitive to magnetic fields with strengths between 0.001 and 0.1 G. 
The reason for this is that the scattering polarization of the 8542 A and 8662 A lines, unlike the 8498 
P^ ' A line, is controlled mainly by the Hanle effect in their (metastable) lower levels. Therefore, in regions 

^Q \ with magnetic strengths sensibly larger than 1 G, their Stokes Q and U profiles are sensitive only to 

the orientation of the magnetic field vector. We also find that for given magnetic field configurations 
the sign of the Q/I and U /I profiles of the 8542 A and 8662 A lines is the same in both atmospheric 
models, while the sign of the linear polarization profile of the 8498 A line turns out to be very sensitive 
O ' to the thermal structure of the lower chromosphere. We suggest that spectropolarimetric observations 

providing information on the relative scattering polarization amplitudes of the Ca ii IR triplet will 
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quiet chromosphere. 
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1. INTRODUCTION 

The chromosphere of the quiet Sun is one of the most complex stehar atmospheric regions (e.g., Harvey 2006, 2009; 
Judge 2006, 2009). Lying between the thin photosphere and the extended 10® K corona, it is here where the magnetic 
field becomes the globally dominating factor ruling the outer solar atmosphere. Measuring the magnetic field vector 
in the solar chromospheric plasma is, however, notoriously difficult, especially outside sunspots and related active 
regions. Quantitative information can be obtained through spectropolarimetry, but the measurement and physical 
interpretation of the weak polarization signatures that the Hanle and Zeeman effects produce in the few chromospheric 
lines that can be observed from ground-based telescopes is not an easy task (see Trujillo Bueno 2010, for a recent 
review) . It is thus important to develop novel methods of "measuring" the chromospheric magnetic field, ideally based 
on the action of the Hanle and Zeeman effects in spectral lines whose line-center intensity images reveal the fine-scale 
thermodynamic structuring of the solar chromosphere. Here we show that the scattering polarization observed in the 
IR triplet of Ca ii contains valuable information on the thermal and magnetic structure of the "quiet" chromosphere, 
from the shock dominated region of the "lower chromosphere" to the "upper chromosphere" below the interface region 
to the 10® K solar corona. 

The only way to obtain direct empirical information on the intensity and topology of the magnetic fields of the 
solar chromosphere is via the measurement and interpretation of polarization signals in chromospheric spectral lines. 
In regions with high concentrations of magnetic flux, such as in sunspots and plages, the polarization patterns are 
dominated by the Zeeman effect. Hence, diagnostic techniques based on this effect are quite useful (e.g., Socas-Navarro, 
Trujillo Bueno & Ruiz Cobo 2000; Socas-Navarro 2005). However, in the quiet Sun, which covers most of the solar disk 
at any given time during the solar magnetic activity cycle, the Zeeman splitting between the tt (AM = M„ — Mi = 0) 
and ab^r (AM = ±1) components (M being the magnetic quantum number) is only a small fraction of the width 
of the spectral lines formed in the chromosphere. As a result, the transverse Zeeman effect makes a very small or 
insignificant contribution to the emergent linear polarization amplitudes. On the other hand, the circular polarization 
of the Zeeman effect as a diagnostic tool is of limited practical interest for the exploration of the magnetism of the upper 
chromosphere because the response function of the emergent Stokes V profiles of strong lines like Ha and Ca ii 8542 
A is significant only in the photosphere and/or lower chromosphere (Socas-Navarro & Uitenbroek 2004; Uitenbroek 
2006). 

There is yet another mechanism producing linear polarization in the spectral lines that can be exploited to diagnose 
the solar chromosphere. Atoms align due to the optical pumping caused by the incident anisotropic radiation field; i.e., 
the individual magnetic M sublevels of energy levels with total angular momentum J > 1/2 are unevenly populated 
(in such a way that the populations of substates with different values of |M| arc unequal) and coherences between them 
may appear. This, in turn, gives rise to a sizable amount of linear polarization in chromospheric spectral lines, either 
by selective emission or absorption of polarization components in the spectral line. A remarkable example is provided 
by the IR triplet of Ca ii. Its linear polarization pattern was considered "enigmatic" for a while (Stenflo, Keller 
& Gandorfer 2000) because it was impossible to achieve even a qualitative fit to the observed relative polarization 
amplitudes considering that scattering line polarization results only from selective emission processes caused by the 
population imbalances of the line's upper level. It was later realized that selective absorption of polarization components 
caused by the presence of populations imbalances in the line's lower level is a very efficient mechanism for producing 
linear polarization in the spectral lines of a stellar atmosphere (Trujillo Bueno & Landi Degl'Innocenti 1997; Trujillo 
Bueno 1999; Trujillo Bueno et al. 2002). In fact, Manso Sainz & Trujillo Bueno (2003a) demonstrated quantitatively 
that the physical origin of the linear polarization pattern observed in the Ca ii lines at 8542 and 8662 A is indeed 
"zero-field" dichroism, i.e., selective absorption of polarization components caused by the presence of atomic alignment 
in the metastable lower levels 4^D3/2 and 4^D5/2 (see Fig. 1), while that observed in the 8498 A line is controled by 
both dichroism and by the selective emission or polarization components that result from the population imbalances 
between the sublevels of its upper level ^P3/2. 

The alignment of the atomic J-levels and the ensuing emergent linear polarization is modified by the presence of 
a magnetic field through the Hanle effect, even if the magnetic field turns out to have a complex, tangled topology, 
with mixed polarities at subresolution scales (e.g., Landi Degl'Innocenti & Landolfi 2004). In the absence of collisional 
depolarization the magnetic field strength B (in G) that is sufficient to produce a significant change is 

^^LlSTxlO-^ (1) 

tmcgj 

where tufc and gj are, respectively, the lifetime (in seconds) and the Lande factor of the J level under consideration. 
Interestingly enough, since the lifetimes of the metastable lower levels of the Ca ii IR triplet are about two orders of 
magnitude larger than the upper level lifetimes and levels with J = 1/2 or J = cannot carry atomic alignment, the 
linear polarization of the Ca ii line at 8662 A (whose upper level has J„ = 1/2) is expected to be sensitive to mC 
fields only, while the scattering polarization of the 8498 A and 8542 A lines (which share the same upper level whose 
Ju = 3/2) can in principle be sensitive also to magnetic fields in the gauss range. In the quiet solar chromosphere 
depolarization of the Ca ii atomic levels by elastic collisions with neutral hydrogen atoms is insignificant for the upper 
levels of the H and K lines, while such collisions have a rather small impact on the atomic polarization of the metastable 
lower levels (Manso Sainz & Trujillo Bueno 2003a; Derouich, Trujillo Bueno & Manso Sainz 2007). Actually, the mere 
fact that the Ca ii 8662 A line shows a clear Q/ 1 polarization signal (see Stenflo et al. 2000) implies that its lower level 



(whose Ji — 3/2) is significantly polarized in the quiet solar chromosphere (Manso Sainz & Trujillo Bueno 2003a). 

The aim of this paper is to investigate the thermal and magnetic sensitivity of the scattering polarization Q / 1 and 
U /I profiles of the Ca ii IR triplet. After formulating in §2 the multilevel problem of scattering polarization and the 
Hanle effect, §4 discusses the unmagnetized case in the "hot" and "cool" chromospheric models described in §3. The 
impact of the Hanle effect on the linear polarization profiles of the Ca ii IR triplet is discussed in §5, showing for each 
solar atmospheric model detailed results for the cases of disk center (§5.1) and close to the limb observations (§5.2). 
The ensuing information is summarized in §5.3 in the form of Hanle effect diagrams for each of the spectral lines of 
the IR triplet. Finally, §6 summarizes our main conclusions with an outlook to future research. 

2. FORMULATION OF THE PROBLEM AND RELEVANT EQUATIONS 

The relevant radiative transfer problem here implies calculating, at each spatial point of the (generally magnetized) 
astrophysical plasma model under consideration, the values of the diagonal and non-diagonal elements of the atomic 
density matrix corresponding to each atomic level of total angular momentum J, which quantify its overall population, 
the population imbalances between its magnetic sublevels, and the quantum coherences between each pair of them. The 
values of such density-matrix elements have to be consistent with the intensity, polarization and symmetry properties 
of the radiation field generated within the medium. This requires solving jointly the RT equations for the Stokes 
parameters and the statistical equilibrium equations for the density-matrix elements. Once such a self-consistent 
solution is obtained it is straightforward to compute the emergent Stokes profiles via the formal solution of the Stokes- 
vector transfer equation for any desired line of sight (LOS). The aim of this section is to formulate such equations with 
emphasis on the case of the IR triplet of Ca ii. 

2.1. The multipolar components of the atomic density matrix 

The excitation state of the atomic system is described through the density matrix (Fano 1957). We consider the 
populations {pj{M, M)), and coherences {pj{M, M')) among the (2 J -I- 1) magnetic sublevels M of each given atomic 
level with total angular momentum J . Alternatively, we use the multipolar components of the atomic density matrix 
(Omont 1977): 

Pq(^)= E(-l)'"''^/2irTT(;^ J^j- ^ )pj{M,m'), (2) 

MM' ^ ^ 

where K = 0, 1, ..., 2 J, Q = —K, —K + 1, ...,K — 1,K, and the expression in brackets is the Wigner 3j-symbol (e.g.. 
Brink & Satchler 1968). The Pq elements with K^O are real numbers measuring population imbalances between 
sublevels, while y/2J + Ip" is the total population of the level. The Pq^q elements are combinations of the complex 
non-diagonal components pj{M, M'). Hermiticity of the density matrix {pj{M,M') — pj{M', M)*) implies p^Q ~ 
{—^)[Pq]* {* meaning complex conjugate). In the absence of net circular polarization in the incident radiation 
Pj{M, M) — pj{—M, —M), i.e., atoms align but they do not orientate. Hence, odd-K pQ elements vanish. Therefore, 
the 5-level atomic model of Ca ii considered here (Fig. 1) requires 1 + 6 + 15 + 1 + 6 = 29 pQ elements to fully 
characterize the excitation state of the atomic system. 

2.2. The transfer equation for the Stokes parameters 
The transfer of spectral line polarization is described by the following equations for the Stokes parameters 

(3b) 
(3c) 

(3d) 

In the Hanle effect regime the magnetic splitting of a and tt components of the line transitions arc negligible and the 
spectral line emissivity thus depends exclusively on the excitation state of the upper level of the transition (i.e., on the 
Pq{Ju) values). In particular, linearly polarized emission is due to the presence of alignment and coherences between 

atomic sublevels of the upper level of the transition. Moreover, since levels are not oriented, ey"° = 0. More explicitly 
(Landi Dcgl'Innocenti & Landolfi 2004; see also Manso Sainz & Trujillo Bueno 1999): 

e'r = eopg + eo^«S!j,{^(3Ai' - l)Po - V2,p^l - p^icosxMpW ~ sinxMpl]) 

+ ^(1 - Ai')(cos2xRcb2] - sin2xlm[p2])}, (4a) 
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l + /i')(cos2xRe[p2]-sin2xIni[p2])|, (4b) 

/!'/,^{v/l-M'(sinxRe[p?] + cosxMpI]) + Ks^n2xRe[pl] + cos2xIni[p2])|, (4c) 

where eg = (hv / 4:^) Aui(f>x-^ V^Ju + 1 (with A/" the total number of atoms per unit volume, Aui the Einstein coefficient 
for spontaneous emission and (j)x the Voigt profile, with x — vq " ^^ and vq the transition's frequency), 6 — arccos (/i) 

(2) 

and X are the inclination and azimuth of the ray, and Wj j is a numerical coefficient that depends only on J„ and 
Ji In particular, w^L ^,^ = 0, w\^L ^..^ — — 2\/2/5, and ur^L ^.^ — \/2/10 (see Table 10.1 in Landi Degl'Innocenti & 
Landolfi 2004). 
The absorption coefficient 777™'', and the dichroism coefficients 77^"° and if^^ are given by analogous expressions but 

with 770 = {hv / A-k)Biu4>x-N' ^f^Ji + 1 instead of eo, Wj j^ instead of Wj^j (wg^ 1/2 — v2/2, u^5/2 3/2 ^ v^/S), and 
with the Pq values of the lower level of the line transition instead of those of the upper level (we neglect stimulated 

emissions). Likewise, since we are assuming that the atomic levels are not oriented, ?7y"° = (when neglecting the 
wavelength shifts between the a components caused by the Zeeman effect). Finally, the anomalous dispersion terms 
(Pq^u) have the same functional dependence on the Pq values of the lower level, but with an antisymmetric dispersive 
dependence on frequency, i.e., Pq^{j = ''lc)''u'4'i^/4'v, with ip^ the dispersion profile. 

Additionally, we consider an unpolarized background continuum providing an extra absorption coefficient 7^™"* and 
emissivity ej°"* = r]f°'^^B^ (where B^ is the Planck function), a very good approximation toward the red part of 
the solar spectrum (e.g., Gandorfer 2000). Therefore, in Eqs. (3), 77/ = 77^"® + 77^™*, e/ = e'j™ + e^™*, while aU 
other emission, dichroism, and anomalous dispersion coefficients have only line contributions (eg — eg"', eu — e}}"", 
TjQ — rfQ°, r]u — 77^"'', pQ — Pq"°, Pu — Pc/"°)- Note that eg and 77g depend on both the population imbalances (pg) 
and on the coherences (pg, with Q = 1,2), while eij and rju depend only on the quantum coherence^ 

The Po elements produce the dominant contribution to the Stokes / profile, while the pg elements (the alignment 
components) contribute to the linear polarization signals, which we quantify by the Stokes parameters Q and U . 
Although ey = 77y = (because, as mentioned above, Pq = 0) according to Eq. (3d) it is possible to generate Stokes 
V in the spectral line due to the —puQ + PqU (second-order) terms. This contribution is however negligible compared 
to the longitudinal Zeeman effect signal produced by the splitting of the cr-components of the emission profile which 
we do not consider here. 

Since the degree of anisotropy of the radiation field in solar-like atmospheres is weak, the population imbalances 
and coherences of each J-level are only a small fraction of its overall population. Therefore, second-order terms can 
be safely neglected in the Stokes- vector transfer equation and the relevant equations for calculating Q/I and U/I are 

-—I « £/ - r]il, (5a) 

as 

as 

^U « [eu - ml] - ViU. (5c) 

as 

Nevertheless, all the calculations of this paper were carried out by solving numerically the full Stokes vector transfer 
equations (3) through the application of the formal solver mentioned in Appendix A. The complication with respect 
to the unmagnetized case considered by Manso Sainz & Trujillo Bueno (2003a) is that now the presence of a magnetic 
field of given strength B, inclination Og and azimuth xb at each grid point of the stellar atmosphere model under 
consideration implies that we have to consider also the quantum coherences (pg, with Qy^O) between each pair of 

magnetic substates pertaining to each J level, in addition to their individual populations (pg)- For this reason, not 
only Stokes Q but also Stokes U is now non-zero, in general. 

The expressions above for the components of the emission vector and propagation matrix (Eqs. 4 and paragraphs 
following them) simplify considerably for several cases of practical interest. For a one-dimensional atmosphere, either 
unmagnetized, or permeated by a microturbulent and isotropically distributed magnetic field of strength _B, or in the 
presence of a deterministic magnetic field with a fixed inclination but with a random azimuth below the line's photon 
mean-free-path, or in the presence of a deterministic magnetic field with a fixed orientation but with a strength in 
the saturation regime of the Hanle effect, the quantum coherences (i.e. the real and imaginary parts of the pf and p^ 
components) vanish in a reference system whose quantization axis must be the parallel to the magnetic field vector for 
the deterministic field case and the the normal to the stellar surface for the three other cases: 

"^ In these expressions, the reference direction for Stokes Q is the perpendicular to the plane formed by the propagation direction and the 
quantization Z-axis of total angular momentum 
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where the p^ {K = 0, 2) values are those of the upper level of the line transition under consideration, and where the 
orientation of the ray is specified by the azimuthal angle x ^'^d by /i = cosO (with 9 the angle between the direction of 
the radiation beam and the quantization axis, which must be the magnetic field direction itself for the deterministic 
field case and the normal to the surface of the stellar atmospheric model for the two other cases). Likewise, the only 
non-zero elements of the propagation matrix are ?]/ and rjQ which are given by identical expressions (i.e. by rji = ej 

and rjQ = eg), but with 770 instead of eo, Wj j instead of Wj j and with the Pq values of the lower level of the line 
transition instead of those of the upper level (for the case in which stimulated emissions are neglected) . 

Under such assumptions, and taking into account the low polarization level in solar- like atmospheres (i.e., that 
Vq/vi ^ 1 ^n.d eg/e/ <^ 1), an approximate formula can be derived to estimate the emergent fractional linear 
polarization amplitude at the center of a strong spectral line (Trujillo Bueno 2003): 
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where Cq = Pg/Po must be evaluated at height in the model atmosphere where the line-center optical distance along the 
line of sight is unity. Its first term is due to selective emission of polarization components (caused by the population 

imbalances of the upper level) , while its second term accounts for "zero-field" dichroism. Note that since w^ 



and w 



(2) 
3/2,1/2 



(2) 
/2,3/2 



= 



V^/2 the scattering polarization of A8662 is exclusively due to "zero-field" dichroism. 



2.3. The statistical equilibrium equations 

General equations for the multivel atom model which neglects quantum interferences between the sublevels pertaining 
to different J- levels are derived in Landi Degl'Innocenti & Landolfi (2004) . For the particular case of the 5- level atom 
of Fig. 1 the following rate equations for the density matrix elements follow: 
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In the previous equations Aut and i?£„ are the Einstein coefficients for spontaneous emission and absorption, re- 
spectively (see Table 1). Jq^ {K^ — 0,2; Q^ — —Kj.,.-.,Kr) are the radiation field tensors introduced by Landi 
Degl'Innocenti (1984). They are integrals over frequency and solid angle of the Stokes parameters (for their explicit 
expressions appropiate for computational work see Manso Sainz & Trujillo Bueno 1999). In particular, Jq is the mean 
intensity integrated over the absorption profile, while 
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The contribution of other Jq components is very small here, where only the magnetic field breaks the symmetry of 
the radiation field. Thus, IJqjq/'^oI — 10^^ throughout the whole outer atmosphere in the five transitions relevant 
for our problem. 

Ciu and Cue are the excitation and desexcitation collisional rates, respectively (see Table 1), between levels i and 

(2) (2) 

u, while C^J and C^/ are collisional transfer rates for alignment between polarizable levels 2, 3, and 5. Derouich et 

al. (2007) have computed values for these transitions in Ca 11. Alternatively one can follow the approach suggested 

(2) 
by Landi Degl'Innocenti & Landolfi (2004), and consider C^^ = CguP where p is just a numerical factor depending on 

the quantum numbers of the transition. D\ is the depolarization rate of the iiT-th multipole of level i due to elastic 
collisions with neutral hydrogen. Theoretical values for Ca 11 are given by Derouich et al. (2007), who confirmed that 

such Dl terms of the statistical equilibrium equations have only a very small impact on the atomic polarization of 



TABLE 1 

Data for 5-level atomic model of Ca ii 



A (A ) u2_ e^. -4„f {s-'-)^ C^i/Ne (s-icm^;[^ 

~ ~ 3000 K 5000 K 7000 K 9000 K 

Allowed transitions 

K 5 i 1.4 X 10** 3.60x10-'^ Z92xkP^ 2.59x10-'^ 2.39x10-'^ 

H 4 1 1.4x10* 3.60x10-^ 2.92x10"^ 2.59x10-'^ 2.39x10-'^ 

8498 5 2 1.11x10*^ 3.68x10"^ 2.93x10"^ 2.51x10-'^ 2.23x10"^ 

8542 5 3 9.6 X 10^* 1.59x10"^ 1.29x10-^ 1.12x10-6 1.01x10"*' 

8662 4 2 1.06x10^ 1.61x10'^ 1.31x10'^ 1.15x10"'^ 1.04x10''^ 

Forbidcn transitions 
~ 5 4 ~ 3.94x10-^ 3.05x10-^ 2.58x10"'^ 2.27x10-* 

— 4 3 — 2.69x10-^ 2.08x10-^ 1.76x10"'^ 1.55x10-'^ 

— 3 1 — 2.20x10"^ 1.71x10-^ 1.44x10-'^ 1.27x10"'^ 

— 3 2 — 5.51x10"^ 4.27x10^^ 3.61x10"'^ 3.19x10"'^ 

— 21 — 2.20x10-^ 1.71x10-^ 1.44x10-'^ 1.27x10"'^ 

''Upper level level of the transition labeled in Fig. 1 

''Lowrer level level of the transition in Fig. 1 

"^Froni NIST atomic spectra database http : //www . nist . gov/physl ab/data/asd ■ cf m | 
"^Froni Shine & Linsky (1974), with A'e being the electron density. 

the metastable lower levels (cf., Manso Sainz & Trujillo Bueno 2003a). 

The first term on the r.h.s. of Eqs. p2|) . (ITill - dTS]) . and (TTSl) represents the relaxation of coherences in the presence 
of a magnetic field with Larnior frequency t'L — 1.3996x10^5 {B expressed in gauss and j^l in Hz), and gj the Lande 
factor of the J-level under consideration. Equations (J10p - (|18p are expressed in a reference system with the quantization 
axis along the magnetic field. The radiative and collisional terms remain formally invariant in any other reference 
system — although clearly, their actual numerical values may change due to the different values of Jq in the new 
reference system. In an arbitrary reference system where the magnetic field has an inclination 9b with respect to the 
quantization axis and an azimuth xb, the magnetic relaxation term reads: 
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where the term multiplying Pg+i {pq_i) is not considered when Q — K {Q = ~K). 

The statistical equilibrium equations (fT0| - (|18p have been derived under several hypotheses. Collisions are considered 
under the impact approximation. The interaction time of the collision is negligible and hence, collisional and radiative 
terms can be added independently. The colliders are assumed to be isotropically distributed (they are Maxwellian) , 
and they cannot generate atomic polarization, though they can transfer it between levelfij. Radiative rates have been 
derived within the framework of the quantum theory of spectral line polarization described in Landi Degl'Innocenti 
and Landolfl (2004), assuming that the radiation field is spectrally flat (i.e., lacking spectral structure), in spectral 
ranges of the order of the energy separation between levels with coherences. This is an excellent approximation for the 
case under consideration since we are considering relatively weak fields (less than ^100 G) and we neglect coherences 
between different J levels. Neglecting J-level interferences is a good approximation for modeling the Ca ii IR triplet 
and the core of the H- and K-lines, but it cannot account for the general Q/ 1 pattern of the UV doublet (see Stenflo 
1980), which is beyond the scope of this work. Finally, we neglect Doppler correlations due to the thermal movement 
of the atoms. All these hypotheses together account for what is called complete frequency redistribution (CRD) in the 
laboraroty frame, since correlations between the incident and scattered radiation fields are completly neglected. 

We calculate the 29 pg unknowns at each grid point of the chosen stellar atmosphere model assuming statistical 
equilibrium 

d 



f^PQ^J) = 0- (21) 

The resulting system of equations is linearly dependent and we substitute one of the equations, say, the one for the 
ground-level population, by the trace equation of the density matrix, which establishes the conservation of the number 
of particles 

^V2J^ + lPo(J«) - 1. (22) 

We have developed a general multilevel radiative transfer computer program for the numerical solution of the non- 
LTE problem of the second kind, which we have applied in this investigation to obtain the self-consistent solution of 
the previous equations. See Appendix A for details. 

^ See Manso Sainz & Trujillo Bueno (2009) for an interesting consequence of this fact for the case of permitted lines at EUV wavelengths 
(i.e., a polarization mechanism of EUV coronal lines and the possibility of mapping the magnetic fields of coronal loops). 



3. THE ATMOSPHERIC AND THE ATOMIC MODELS 

High spatial and temporal resolution intensity images taken at the line center of strong absorption lines, like Ha 
and Ca ii 8542 A, show that the chromosphere of the quiet Sun is a highly inhomogeneous and dynamic medium (e.g., 
review by Rutten 2007). The spatial and dynamic complexity is such that it is obvious that the two one-dimensional 
(ID) atmospheric models we have chosen for this investigation (i.e., the "hot" FAL-C model of Fontenla et al. 1993 
and the "cool" M-CO model of Avrett 1995) should only be considered as rough representations of the stratification of 
the kinetic temperature and density at two different phases of inter-network oscillations (see Fig. 2). Nevertheless, as 
we shall see below, our solution of the multilevel Hanle-efFect problem for the IR triplet of Ca ii in such semi-empirical 
models of the solar atmosphere allows us to reach the main goal of this paper, namely to demonstrate the diagnostic 
potential of the linear polarization signals produced by atomic polarization in the levels of the Ca ii IR triplet. 

We consider a 5-level atomic model of Ca ii including the ground level, the metastable ^1)3/2,5/2 levels and the 
^Pi/2,3/2 upper levels of the H and K lines (see Fig. 1). The excitation state in such a model is given by 29 statistical 
tensor elements pn (K = 0, 2, 4; Q — ~K, ...,K), including the total populations of the five levels. The radiative 

transition probabilities are those of Edlen & Risberg (1956) compiled at the NIST databascQ. (Note that these 
transition probabilities for the infrared triplet are slightly larger than those used by Manso Sainz & Trujillo Bueno 
(2003a) for the unmagnetized case; correspondingly, the Q/I values for the zero-field case are slightly larger here than 
in our previous work.) Inelastic collisional rates and line broadening parameters have been taken from Shine & Linsky 
(1974; see Table 1). Elastic depolarizing collisions are treated according to Lamb & ter Haar (1971), which gives 
results similar to those obtained using the collisional rates given by Derouich et al. (2007). The collisional pumping 
oi K y^ statistical tensors are treated following Landi Degl'Innocenti & Landolfi (2004). 

We have calculated the number density of Ca 11 ions at each atmospheric height of the model under consideration by 
solving the standard non-LTE radiative transfer problem for a realistic atomic model that includes also the bound-free 
transitions from all the bound levels of Ca 11 to the ground level of Ca III. With the resulting number densities of Ca 
II ions we then solved the non-LTE problem of the second kind for the atomic model of Fig. 1 (see the Appendix), 
which is sufficiently realistic for investigating the polarization signatures produced by atomic level polarization and 
the Hanle effect in the lines of the Ca 11 IR triplet. We point out that while Bh (see Eq. 1) is of the order of a few 
gauss for the upper level of the 8498 A and 8542 A lines, it lies in the milligauss range for the lower (metastable) levels 
of the Ca 11 IR triplet. 

4. THE UNMAGNETIZED REFERENCE CASE 

For symmetry reasons, in a one-dimensional model atmosphere without magnetic fields the only non-zero Stokes 
parameters are /(A) and Q{X) (choosing the polarization directions as in Eqs. (jlj). 

Figure 3 shows the wavelength variation of the emergent fractional linear polarization {Q/I) in the Ca 11 IR triplet, 
calculated for a line of sight with /x — 0.1 in the "hot" FAL-C model (middle panels) and in the "cool" M-CO model 
(bottom panels). While Q/I in the 8542 A and 8662 A lines are positive in both models, the Q/I profile of the 8498 
A line is positive in the "hot" model but negative in the "cool" model. Moreover, the 8498 A Q/I profiles calculated 
in such atmospheric models have even different shapes, and their amplitudes are one order of magnitude smaller than 
those of the 8542 A and 8662 A lines. The (positive) Q/I amplitudes of the 8542 A and 8662 A lines turn out to be 
significantly larger in the "cool" model. This dependence of the amplitude of the Q/I profiles on the chosen model 
atmosphere is not at all surprising, because the atomic polarization that anisotropic pumping processes induce in the 
atomic energy levels sensitively depends on the temperature structure of the model. In particular, the steeper the 
temperature gradient, the larger the anisotropy factor and hence, the induced atomic level polarization (e.g., see figure 
4 in Trujillo Bueno 2001). The upper panels of Fig. 4 show the spatial variation of the "degree of anisotropy" in each 
of the considered atmospheric models, indicating for each LOS characterized by p, the atmospheric height where the 
line-center optical depth is unity along the LOS. Note that at the height corresponding to a LOS with p, = 0.1 the 
anisotropy is larger in the M-CO model. The decrease of the anisotropy factor at very low optical depths is due to 
the rapid widening of the (px profile with which the radiation field is averaged (e.g., see Eq. (18) for the definition of 
the Jq radiation field tensor component). This is a consequence of the temperature rising in the upper chromosphere- 
transition region. The two lower panels of Fig. 4 show the behavior of the fractional atomic alignment of the levels 
of Fig. 1 having J > 1/2 (i.e., the two metastable levels, which are the lower levels of the Ca 11 IR triplet, and the 
upper level of the K-line). The fact that the sign of the Q/I profile of A8498 (i.e., the weakest line of the Ca 11 IR 
triplet) is very sensitive to the thermal structure of the lower chromosphere can be easily understood by combining 
the information provided in the lower panels of Fig. 4 and the approximate Eddington-Barbier formula (8). 

Fig. 3 suggests that the scattering polarization of the Ca 11 IR triplet can be used as a sensitive thermometer of the 
"quiet" regions of the solar chromosphere, as already noted by Manso Sainz & Trujillo Bueno (2001). In this respect, 
a particularly useful quantity is the center-limb variation of the Q/I line-center amplitudes (see Fig. 5). As seen in 
this figure, the most noteworthy feature is that in the "hot" FAL-C model the (positive) Q/I amplitude of the 8498 A 
line is larger at /i«0.2 than at /i«0.1 and that it becomes negative for line of sights with /i values significantly smaller 
than 0.1. 

5. THE HANLE EFFECT IN THE CA II IR TRIPLET 

^ |http : //physics ■ nlst . gov/PhysRef Data/ ASP | 



The aim of this section is to show how the Hanle effect modifies the emergent Unear polarization of the IR triplet of 
Ca II. To this end, we consider the following two geometries. 

5.1. Forward scattering geometry 

In the absence of magnetic fields and horizontal atmospheric inhomogeneities the polarization of the atomic levels do 
not produce any linear polarization in the spectral line radiation observed at disc center (/i = 1; see Fig. 5). The same 
is true in the presence of a vertical magnetic field, because in a one-dimensional stellar atmospheric model the vertical 
direction coincides with the symmetry axis of the anisotropic radiation field that induces atomic level polarization. 
However, in the presence of an inclined magnetic field the symmetry of the problem is broken and, as a result, the 
ensuing atomic level polarization can generate linear polarization even for a line of sight with /i = 1. In this case, the 
linear polarization is created by the Hanle effect of the inclined field (see a demonstrative observational example in 
Trujillo Bueno et al. 2002). 

Figure 6 shows the fractional linear polarization of the IR triplet of Ca ii generated in forward scattering geometry 
by the Hanle effect of a horizontal magnetic field in the FAL-C (top panels) and M-CO (bottom panels) models. In 
the ("hot") FAL-C model, the linear polarization amplitudes are at the level of 10""* when the magnetic strength of 
the horizontal field is similar or larger than only 0.1 gauss, with Q/I > (i.e., parallel to the horizontal magnetic 
field) for the 8498 A line and with Q/I < (i.e., perpendicular to the horizontal magnetic field) for the 8542 A and 
8662 A lines. In the ("cool") M-CO model Q/I > in the 8542 A and 8662 A lines, with amplitudes again at the 
10~^ level for _B>0.1 gauss, while the shape of Q/I in the 8498 A line is very peculiar (i.e., similar to the typical 
observational signature of the transverse Zeeman effect) and with its maximum amplitude one order of magnitude 
smaller. Figure 7 shows the results of our RT calculations for the case of a magnetic field inclined by only 30°. Note 
that now, U/Iy^O. At present, detection of these weak linear polarization signals is challenging, though feasible. But 
they will become straigtforward with future large-aperture solar telescopes such as the Advanced Technology Solar 
Telescope (Keil, Rimmele & Wagner 2009) or the European Solar Telescope (CoUados 2008). 

5.2. Close to the limb geometry 

In forward scattering geometry, symmetry imposes that linear polarization vanishes if the magnetic field is inclined 
but has a random azimuth distribution, unlike in the fixed inclination and azimuth case just considered. However, for 
observations away from the solar disk center (/x<l), scattering may produce linear polarization in the emergent spectral 
line radiation even for magnetic fields with a uniformly distributed azimuth within the spatio-temporal resolution 
element of the observations. This is illustrated in Fig. 8, which shows the linear polarization amplitudes of the Ca ii 
IR triplet in the FAL-C model (top panels) and in the M-CO model (bottom panels) vary with the magnetic strength 
(B) and the inclination {0b) of the assumed random-azimuth field for a close-to-the-limb observation (/x = 0.1). 

Figure 8 shows how sensitive the Q/I line-center amplitude of the 8498 A line is to the model's thermal structure, 
and its very peculiar behavior in the sub-gauss range. The 8542 A and 8662 A lines are mostly sensitive to magnetic 
fields between 0.001 and 0.1 gauss, although the 8542 A line is also weakly sensitive to fields between 1 and 10 gauss. 
The 8542 A line enters into the saturation regime of the upper-level Hanle effect for fields stronger than ^10 gauss, 
while just 0.1 gauss saturate the 8662 A line. On the contrary, the 8498 A line is highly sensitive in the whole regime 
0.001-10 gauss. As seen in Fig. 3, in the absence of magnetic fields the linear polarization amplitude of the 8498 A 
line for a LOS with /i = 0.1 is one order of magnitude smaller than those corresponding to the 8542 A and 8662 A 
lines. However, the three Ca ii lines may show similar Q/I amplitudes if there is a magnetic field of the order of only 
0.1 gauss in the atmospheric region where the linear polarization of the 8498 A line originates. 

Figure 9 shows the emergent Q/I amplitudes averaging the Q/I profiles of both atmospheric models with a weight 
of 0.55 for the FAL-C profiles and 0.45 for the M-CO profiles. These weights are similar to those used by Avrett (1995) 
to simultaneously reproduce the intensity profiles observed in the 4.6 /im lines of CO and in the H & K lines of Ca II, 
and to those used by Holzreuter et al. (2006) to fit the Q/I profile of the Ca ii K-line observed in a quiet regions very 
close to the solar limb. These type of fits obtained by mixing a cool and a hot component should not be considered 
as a realistic representation of the chromospheric conditions, but as further indication of the complexity of the (time- 
dependent) three-dimensional thermal structure of the solar chromosphere. Nevertheless, this two-component model 
shows that the Q/I observations of the IR triplet of Ca ii reported by Stenflo et al. (2000) are compatible with very 
weak magnetic fields in the "region of formation" of the 8498 A line (i.e., in the lower chromosphere) and with fields 
stronger than 10 gauss (and inclined by about 30°) in the "regions of formation" of the 8542 A and 8662 A lines (i.e., 
in the upper chromosphere). 

Stokes U is zero if the magnetic field has a random-azimuth distribution, and it is non-zero if the atmosphere is 
permeated by a magnetic field vector with a well-defined inclination and azimuth. This is illustrated in the upper 
panels of Fig. 10, which corresponds to the FAL-C model. In the absence of magnetic fields, only Q/I is non-zero (see 
the dashed lines), but U/I becomes significant even at fields strengths as low as 0.005 gauss. Note that for this case 
(B=0.005 gauss), characterized by a horizontal magnetic field pointing towards the observer (triangles), perpendicular 
to the LOS (squares) and away from the LOS (diamonds), the Stokes Q/I profiles of the Ca ii IR triplet lines are 
always positive, while U/I changes its sign. The dotted lines of Fig. 10 show the case of a 100 gauss horizontal field 
pointing towards the observer (i.e., a strength for which the three IR lines of Ca ii are in the saturation regime of the 
upper-level Hanle effect). The bottom panels of Fig. 10 show the results of similar calculations for the cool M-CO 
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model. In this model the Q / 1 profile of the 8498 A line changes its sign and shape when the magnetic field azimuth 
changes from xb ~ 0° (or xb = 180°) to xb ~ 90° 

5.3. Hanle-effect diagrams 

A Hanle-effect diagram shows the line-center amplitudes of Q/I and U/I varying with the strength and the azimuth 
of the magnetic field corresponding to a given line-of-sight and inclination of the magnetic field vector. Figure 11 shows 
restricted Hanle-effect diagrams corresponding to the configurations considered in Fig. 10 (close-to-the-limb /i = 0.1 
observation). In particular, the thin solid lines show how linear polarization varies with the magnetic strength, for a 
horizontal magnetic field pointing towards the observer. We recall that 50 gauss are needed to completely depolarize 
the 8498 A line, while just 10 gauss are sufficient for the 8542 A line and only 0.1 gauss for the 8662 A line (see Fig. 
8). Thick solid lines show, for the case of a 0.005 gauss horizontal field, how the Q/I and U/I line center amplitudes 
vary as the orientation of the magnetic field vector changes from pointing towads the observer (triangles), to being 
perpendicular to the line of sight (squares) to pointing away from the observer (diamonds). Comparing the upper and 
lower panels we see the differences due to the different thermal structures of the FAL-C and M-CO models. 

Figures 12 and 13 show full Hanle-effect diagrams obtained from the emergent Stokes Q/I and U/I profiles calculated 
for a close-to-the-limb (/i = 0.1) line of sight in the "hot" and "cool" atmospheric models, respectively. Each panel 
corresponds to a fixed inclination of the magnetic field vector. Each of the four thick lines in any given panel 
correspond to a fixed strength of the magnetic field with azimuth ranging from xb = 0° to xb = 180° (solid line), 
and from XB = 180° to xb = 360° (dashed line). Thin lines correspond to magnetic fields ranging between G and 
1000 G for a fixed xb- They show the pattern of depolarization {Q/I decreases) and rotation of the direction of linear 
polarization {U/I appears) characteristic of the Hanle effect in a close to the limb scattering geometry. 

Finally, Figs. 14 and 15 show Hanle-effect diagrams for a disc center observation (/x = 1). Note that in this forward 
scattering geometry the Hanle effect of a magnetic field inclined with respect to the solar local vertical creates linear 
polarization. 

6. CONCLUSIONS 

We have investigated theoretically the sensitivity of the scattering polarization Q/I and U/I profiles of the Ca 
II IR triplet to the presence of magnetic fields through the Hanle effect. To this end, we have applied a multilevel 
radiative transfer code (see Appendix A) for treating the transfer of spectral line polarization due to scattering in 
weakly magnetized stellar atmospheres. We find that the Ca ii 8542 A and 8662 A lines are mostly sensitive to fields 
between 0.001 and 0.1 G, the former being also weakly sensitive to fields between 1 and 10 G. Otherwise the linear 
polarization signals are dependent only on the magnetic field geometry, but not to its intensity. The 8498 A line is 
highly sensitive in the whole regime 0.001-10 G, but its scattering polarization amplitude is significantly lower in the 
absence of magnetic fields. 

Probably, the most interesting line of the Ca ii IR triplet because of its diagnostic potential is the strongest one: 
A8542. Its linear polarization is sensitive to the orientation of the magnetic field in the chromosphere of the quiet Sun, 
though not too much to its strength unless i? < 1 G there. Although its scattering polarization amplitude depends on 
the thermal structure of the solar chromosphere, the emergent Q/I and U/I profiles do not change sign between the 
("hot") FAL-C model and the ("cool") M-CO model, unlike A8498, which is very sensitive to the physical conditions of 
the lower chromosphere. The thermal sensitivity of the A8662 line (which is sensitive only to milligauss field strengths) 
is similar to that of the 8542 A line. Therefore, spatial fiuctuations in the sign of the emergent U/I profiles of the 
Ca II 8542 A and 8662 A lines may be ascribed more safely to spatial variations in the azimuth of the chromosphcric 
magnetic field than is the case for A8498. 

More generally, the emergent Q/I and U/I profiles are produced by the joint action of atomic level polarization, and 
the Hanle and transverse Zeenian effects. In this paper we have considered only the action of atomic level polarization 
and the Hanle effect. These effects dominate the emergent linear polarization profiles for inclined magnetic fields with 
strengths weaker than Bq, where the Bq value depends on the scattering geometry. 

In the forward scattering geometry of a disk center observation (/i = 1) the linear polarization of the Ca ii IR triplet 
is dominated by the Hanle effect if the magnetic field is weaker than about 10 G. In fact, while the contribution of the 
transverse Zeeman effect is negligible for 0<i3<10 G the Hanle effect creates fractional linear polarization signals of 
the order of 10""* already for horizontal magnetic fields as weak as 0.1 G (see Figs. 6 and 7). For magnetic strengths 
10<i3<100 G the contribution of the transverse Zeeman effect to the linear polarization observed at the solar disk 
center should not be neglected. Detection of Q/I and U/I disk center signals caused either by the Hanle effect alone 
(if B<10 G) or by the joint action of the Hanle and transverse Zeeman effects (if 10<i3<100 G) would require very 
high polarimetric sensitivity together with a spatial and temporal resolution sufficient to at least resolve partially the 
magnetic field azimuth. 

Given the weakness of the theoretical scattering polarization signals (^lO""* at /i = 1; <10^'^ at /x = 0.1), the 
most favourable geometry for observing them is in quiet regions close to the limb. In this geometry the emergent 
linear polarization of the Ca ii IR triplet is non-zero even if the magnetic field has a random azimuth within the 
spatio-temporal resolution element of the observation (see Figs. 8 and 9). Moreover, for magnetic strengths sensibly 
weaker than 100 G the contribution of the transverse Zeeman effect to the emergent linear polarization is expected to 
be smaller than that caused by atomic level polarization. For stronger fields (e.g., 100<i?<500 G) the wings of the 
emergent Q/I and U/I profiles may show sizable features produced by the transverse Zeeman effect, but the line-center 
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amplitudes would still be dominated by atomic level polarization. 

Finally, we point out that information on the relative amplitudes and signs of the scattering polarization profiles 
of the Ca ii IR triplet would be very useful to help constrain the thermal and magnetic structure of the quiet solar 
chromosphere. In principle, to determine the structure of the magnetic field via Stokes inversion techniques is possible, 
although it would be easier if one is able to infer first the atmospheric thermal structure from the observed intensity 
spectrum. The Q/I and U /I signals we have studied here should also be exploited to investigate the reliability of 
three-dimensional models of the solar chromosphere resulting from holistic magnetohydrodynamic (MHD) simulations 
of the solar atmosphere. To this end, it suffices with confronting spectropolarimetric observations of the Ca ii IR 
triplet with synthetic Stokes profiles obtained by solving the non-LTE problem of the 2nd kind in 3D snaphot models 
resulting from such MHD simulations. 

APPENDIX 

The Numerical Solution of the non-LTE Problem of the Second Kind 

Modeling the spectral line polarization produced by atomic level polarization and its modification by the Hanle effect 
requires calculating, for multilevel systems, the excitation and ionization state of chemical species of given abundance 
that is consistent with both the intensity and polarization of the radiation field generated within the (generally 
magnetized) plasma under consideration. This so-called non-LTE problem of the 2nd kind (cf., Landi Degl'Innocenti 
& Landolfi 2004) is a very involved non-local and non-linear radiative transfer problem which requires solving the 
statistical equilibrium equations for the elements of the atomic density matrix and the Stokes- vector transfer equation 
for each of the allowed transitions in the multilevel model. Once such a self-consistent excitation state is known 
throughout the medium, it is then straightforward to solve the transfer equation for any desired line of sight in order 
to obtain the emergent Stokes profiles to be compared with spectropolarimetric observations. To this end, we developed 
a general computer program, which we will describe in greater detail in a future publication along with the efficient 
and accurate radiative transfer methods on which it is based (for a first advance see the workshop contribution by 
Manso Sainz & Trujillo Bueno 2003b; see also Manso Sainz 2002). In this Appendix we present only a brief summary 
of the code with emphasis on its application to solve the problem of scattering polarization and the Hanle effect in the 
IR triplet of ionized calcium. 

A summary of the numerical approach is as follows: 

1. The equations are formulated within the spherical tensors representation of the density matrix and radiation 
field tensor. There are (2 J -I- 1){J -\- 1) unknowns Pq with K even, for each level with integer angular momentum 
J, and (2 J + 1) J unknowns for each level with semi-integer angular momentum J. 

2. The radiative transfer equations for the Stokes parameters are integrated along short characteristics using the 
quasi-parabohc DELO method (DELOPAR; Trujillo Bueno 2003). 

3. The iterative corrections for the unknowns (the statistical tensors pg ) are calculated applying a suitable gener- 
alization to the multilevel atom case of the fast iterative methods described in Trujillo Bueno & Manso Sainz 
(1999). This implies writing down the statistical equilibrium equations taking explicitly into account the contri- 
bution of the diagonal components of the A operator and linearizing according to Eqs. (53) and (54) of Trujillo 
Bueno (2003). 

4. The iterative scheme can speeded up further through polynomial Ng-accclcration (e.g., Aucr 1987). 

The number density of the ion under consideration is computed in non-LTE. At each point in the atmosphere we solve 
the statistical equilibrium equations for the multipolar components of the atomic density matrix plus the conservation 
of particles equation to calculate the excitation state of the atomic system. We calculate the radiation field at each 
point in the model atmosphere for each of the allowed radiative transitions in the model atom by formal integration 
of the radiative transfer equations for the Stokes parameters applying the DELOPAR method. Then, radiation field 
tensors Jq are calculated and used to work out the new excitation state of the atomic system and so on, iteratively. 
Figure 16 shows the convergence behavior of the above-mentioned iterative method for the numerical solution of the 
(Ca ii) non-LTE problem of the second kind. 

Financial support by the Spanish Ministry of Science and Innovation through project AYA2007-63881 (Solar Mag- 
netism and High-Precision Spectropolarimetry) is gratefully acknowledged. 
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Fig. 1. — The atomic model for Ca II. Labels indicate the total angular momentum of the level and its numbering (between parenthesis) 
in Eqs. I|10 |l - ljl8p . Solid lines connecting levels show allowed transitions and their labes their wavelength in A. 
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Fig. 2. — The atmospheric models. Variation with height of the kinetic temperature (left panel), hydrogen and electron density (right 
panel) in the FAL-C model (solid lines), and M-CO model (dashed lines). Note in the left panel that the "lower chromosphere" of the 
M-CO model is relatively cool while it is relatively hot in FAL-C. 
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Fig. 3. — The unmagnctized case. Emergent Q/I profiles of the Ca II IR triplet calculated for a line of sight with jj, = 0.1 in the "hot" 
FAL-C model (middle panels) and in the "cool" M-CO model (bottom panels). The positive Q-direction is parallel to the nearest limb. 
Note that the sign of the Q/I profile of the 8498 A line is positive in the "hot" model and negative in the "cool" model. The uppermost 
panel shows the emergent intensity profiles computed in the FAL-C model (solid lines) and in the M-CO model (dotted lines). 
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Fig. 4. — The variation with the (vertical) K-hne optical thickness of the degree of anisotropy in each of the lines of the Ca II model of 
Fig. 1 after obtaining the self-consistent solution in the FAL-C model (left panel) and in the M-CO model (right panel). For each of the 
IR triplet lines the corresponding thick horizontal segment indicate where its line-center optical depth is unity for line of sights going from 
fi = 1 (the lowest tk of each segment) to ^ = 0.1 (the largest tk of each segment). 
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Fig. 5. — The unmagnctizcd case. Ccnter-to-limb variation of the hne-center scattering polarization amplitudes in the Ca II IR triplet 
produced by atomic level polarization in the "hot" FAL-C model (top panels) and in the "cool" M-CO model (bottom panels). The positive 
Q-direction is parallel to the nearest limb. Note that for relatively small ^ values the Qjl of A8498 is negative for both model atmospheres. 
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Fig. 6. — The Hanle effect of a horizontal magnetic field in forward scattering. The emergent Qjl profiles of the Ca II IR triplet calculated 
for a line of sight with /.t = 1 (i.e., disk center observation) assuming the presence of a horizontal magnetic field of 100 G (solid lines), 0.1 G 
(long-dashed line), 0.01 G (short-dashed line), and 0.001 G (dotted line) in the "hot" FAL-C model (top panels) and in the "cool" M-CO 
model (bottom panels). The positive Q-direction is along the magnetic field. 
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Fig. 7. — The Hanle effect of a non- horizontal magnetic field in forward scattering. As in Fig. 6, but for a magnetic field inclined by 30° 
Note that in this case Stokes \J is not zero. 
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Fig. 8. — The Hanlc cfEect for a close to the hmb hne of sight. The emergent Q/I hne-center amphtudes of the Ca II IR triplet calculated 
for a line of sight with /i = 0.1 in the "hot" FAL-C model (top panels) and in the "cool" M-CO model (bottom panels), assuming a 
magnetic field with the indicated inclination and a uniformly distributed azimuth within the spatio-temporal resolution element. The 
positive Q-direction is the parallel to the nearest limb. 
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Fig. 9. — The Hanle effect for a close to the hmb hne of sight. The emergent Qjl hne-ccnter amphtudes of the Ca II IR triplet calculated 
for a line of sight with /.j = 0.1 assuming that 55% of the radiation originates in the "hot" FAL-C model and 45% in the "cool" M-CO model, 
assuming a magnetic field with the indicated inclination and a uniformly distributed azimuth. The positive Q-direction is the parallel to 
the nearest limb. 
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Fig. 10. — The emergent Qjl and V jl profiles of the Ca II IR triplet calculated for a line of sight with /i = 0.1 in the "hot" FAL-C 
model (top panels) and in the "cool" M-CO model (bottom panels), assuming a magnetic field with the following magnetic configurations: 
pure scattering case (dashed lines) and the case of a magnetic field pointing to the observer with B=0.005 G (solid lines with the triangle 
symbol) and B=100 G (dotted lines). Actually, for the B=0.005 G case three different orientations of the linc-of-sight are shown; along 
the field (triangles), perpendicularly to the field (squares), and away from the field (diamonds). The positive Q-dircction is the parallel to 
the nearest limb. 
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Fig. 11. — The evolution of the line center signal in the Q/I - U/I space, for the FAL-C model (upper panels) and the M-CO model 
(bottom panels). Thin lines: horizontal magnetic field with variable strength (from B=0 G to B— >oo). Thick lines: B=0.005 G horizontal 
field observed from diffecrent angles. The symbols correspond to the particular configurations shown in Fig. 10. 
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Fig. 12. — Hanle diagrams showing the Qjl and V jl signal at the line center of the Ca II infrared triplet lines for several line-of-sight 
and magnetic field configurations in the FAL-C atmospheric model. The observation is always close-to-the-limb (^ = 0.1). Thick lines 
correspond to constant strength magnetic fields of 1.8 X 10^^, 5.6 X 10~^, 1.8 X 10~^, and 18 G when observed with 0° < X — Xb < 180° 
(solid line) and 180° <% — XB < 360° (dotted line), where x ~ XB is the azimuth between the line-of-sight and the magnetic field. (These 
two tracks do not coincide in the Qb = 90° case because we are not in the 90° -scattering or tangential observation limit, but at ^ = 0.1). 
The thin solid lines correspond to x — XB = 0°, 30°, 60°, 90°, 120°, 150°, and 180°. Upper panels correspond to a magnetic field inclined 
0s = 90° with respect to the vertical (i.e., horizontal magnetic field), central and lower panels correspond to inclinations of 0g = 60° and 
ds = 30°, respectively. 
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Fig. 13. — Like in Fig. 12 but for the M-CO model atmosphere. 
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Fig. 14. — Like in Fig. 12 but in forward scattering geometry. 
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Fig. 15. — Like in Fig. 13 but in forward scattering geometry. 



0.02 



0.03 L 
0.03 -0.03 



0.02 



0.01 0.00 0.01 

U/l (%) 



0.02 0.03 



27 








20 40 

Iteration Number 



60 



80 



Fig. 16. — Maximum relative change (i?c(it) = Maxi=i_..._jv^ [IPq (*; it) ~ Pni^'i '* ~ 1)I/IPq (*; it)l]) fo'' each of the 29 Pq components in 
a radiative transfer calculation using a 5-level atomic model of Ca II. Solid lines show the convergence rate of all pp components. Heavy 
solid lines correspond to pfj. Note that due to superpositions only three lines are visible while there are 5 pjj components. Dashed lines 
correspond to the case in which Ng-acceleration is used. For clarity, in this case only the pg convergence rates are shown. 



